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Abstract 

We study the random held p-spin model with Ising spins on a fully connected 
graph using the theory of large deviations in this paper. This is a good 
model to study the effect of quenched random held on systems which have a 
sharp hrst order transition in the pure state. For p = 2, the phase-diagram 
of the model, for bimodal distribution of the random held, has been well 
studied and is known to undergo a continuous transition for lower values of 
the random held {h) and a hrst order transition beyond a threshold, htp{^ 
0.439). We hnd the phase diagram of the model, for all p> 2, with bimodal 
random held distribution, using large deviation techniques. We also look at 
the huctuations in the system by calculating the magnetic susceptibility. For 
p = 2, beyond the tricritical point in the regime of hrst order transition, 
we hnd that for htp < h < 0.447, magnetic susceptibility increases rapidly 
(even though it never diverges) as one approaches the transition from the 
high temperature side. On the other hand, for 0.447 < h < 0.5, the high 
temperature behaviour is well described by the Curie-Weiss law. For all 
p > 2, we hnd that for larger magnitudes of the random held {h > ho = 1/p!); 
the system does not show ferromagnetic order even at zero temperature. 
We hnd that the magnetic susceptibility for p > 3 is discontinuous at the 
transition point for h < ho- 
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1. Introduction 


Magnetic systems in the presence of random field are an important class 
of disordered systems which have been extensively studied for last many 
years. These systems are governed by two energies: the exchange interaction 
between spins which want spins to be correlated and local random held h 
at each site which would make the spins uncorrelated. Thermal frustrations 
play a secondary role in random held models, since their disordering ehect is 
much weaker [l|, . 

One of the very well studied model to understand the ehect of random 
held is the random held Ising model (RFIM), which was introduced by Imry 
and Ma in 1975 j^. The pure model (Ising model) has a continuous transition 
for all d > 2. The model has been very useful as it has found its physical 
realizations in many of the condensed matter systems j^. Hence the model 
has attracted a lot of attention, both theoretically and experimentally BH 
S B 0 B B [loj |. The model has many interesting and intriguing properties. 
For example, at the mean held level the phase-diagram of RFIM is known to 
change its character as a function of the random held distribution function. 
For gaussian distribution, the transition stays continuous for all strengths of 
the random held, but for bimodal distribution it exhibits a tricritical point 
0- In three dimensions though the existence of tricritical point has been 
widely debated ji, IS; 4, 13, 13 1, with recent numerical simulations suggesting 
continuous transition even in the strong disorder regimejs], 13 . 


Related model is a p-spin model, in which p-spins are connected by an 


exchange interaction [1^. Pure system has a sharp hrst order transition for 


p > 3, which can easily be studied in the mean held limit. The random held 
version of this model, for p > 3, has relatively attracted much less attention 
in the literature. This model is a good system to study the ehect of local 
held randomness on sharp hrst order transitions [iB- The model is also of 
relevance since it is known to have connections with the structural glasses 
00 . 

Ehect of quenched local randomness for systems which undergo continu¬ 
ous phase transition in the absence of randomness is well described by the 
Harris criterion [13]. This is possible due to the universal behaviour of the 
correlation length. On the other hand, even though hrst order transitions are 
common in nature, ehect of random local disorder on hrst order transitions 
is much harder to study theoretically, as the correlation length remains hnite 


and hence the universal behaviour is not expected in these systems [19 
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In this paper we look at the random held p-spin model on fnlly connected 
graphs using the theory of large deviations. Large deviations lend itself nat¬ 
urally to the study of mean held statistical mechanics models . Recently, 
Lowe et al derived the rate function (to be dehned in section 2) for the 
magnetisation of random held Ising model on fully connected graphs using 
a method based on tilted large deviation techniques j^. Their method is 
much simpler than the earlier approaches and can easily be generalised 
to many other random held problems. In this paper, we extend it to solve 
the random held p-spin model for arbitrary p. Note that even though cal¬ 
culations are a straightforward generalisation of j^, they are useful as they 
allow us to study the random p-spin model fully and hence allow us to look 
into the issue of understanding the disorder ehects near the hrst order tran¬ 
sitions. Besides obtaining the phase diagram in this paper, we also derive 
a method to obtain the value of quantities like magnetic susceptibility from 
the rate function. Magnetic susceptibility is an experimentally measurable 
quantity, and hence is useful to make direct connection with the experimen¬ 
tal observations. We hnd that in the hrst order regime of RFIM, magnetic 
susceptibility is a rapidly increasing function of the temperature as one ap¬ 
proaches the critical point, even though it does not diverge at the transition 
point. For p > 3, it is discontinuous at the transition point, with a typical 
paramagnetic behaviour in the high temperature phase. 

The plan of the paper is as follows: In section HI we give the details of 
the model and calculate the rate function. We also derive the expression 
for magnetic susceptibility in this section. In section |3l we give the phase- 
diagram obtained by looking at the value of magnetisation, and in section |H 
we calculate the magnetic susceptibility for p = 2, 3 and 4. We conclude in 
section O 

2. Rate function for a random field p-spin model 

We consider a Hamiltonian of the form: 

^ N k=p N N 

^ ^ ~ p\NP-^ ^ n X] (^) 

where are Ising spins which take values ±1, H is the uniform external 
field and hi is the site dependent random field taken from a distribution. The 
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probability of getting a particular configuration for a given realization 

of {hi} is given by: 


PN{/3,{siJ,{hi}) =-j^exp I — Y[si^+ /3Yif^ii + H)i 


( 2 ) 


where, 


{hi}) = ^ exp ( ^ Y 


N k=p 


N 




2N 


'll 


( 3 ) 


2i ,....^ip k —1 


21 = 1 


is the partition function of the model for a given realization {hi} of the 
random held and /3 is the inverse of the temperature. In order to get the 
free energy, one needs to average the log of the partition function over the 
distribution of random helds. This averaging over the log is a hard problem, 
and has been done only for a relatively fewer models like the random energy 
model . Another approach is the replica approach where the < log Z > 
is approximated by the use of the identity < Z” >= 1 + n < log Z > +0{n^) 

. One more method which has been found useful to study the mean held 
spin systems is the large deviation technique. One can prove the existence of 
large deviation principle with an explicit rate function for the order parameter 
of the model and then use it to study the phase diagram j^. Recently, for 
random held Ising model (p = 2), Lowe et al derived the rate function 
explictly for the probability distribution of magnetisation. For a sequence of 
random variables {sj}, the rate function is dehned as 


lim e A) = -/(A), (4) 

N^OO iV 

where P/v({sj} G A) is the probability that the given sequence of random 
variables ({sj}), belong to the set A. For our model, A can be dehned via 
X = he., all sequences with same value of x will have the same 

probability, as given by the equation above. 

Note that trying to get the rate function directly, using say Gartner Ellis 
theorem j^, is quite complicated asY^Si/N is not a sum of N independent 
random variables with respect to the probability measure in Eq. |2l In order 
to calculate the rate function Lowe et al developed a method which 
mainly involved the following steps: 
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1. Prove large deviation principle (LDP) with respect to a measure Q, 
which is a measure on non-interacting spins and hence is easier to average 
over the quenched held. 

2. Use titled LDP to calculate the rate function for the probability dis¬ 
tribution of magnetisation in the presence of coupling term from the rate 
function of Q. 

This method turns out to be very useful especially in the cases where 
the randomness is not in the the coupling term. Also it gives the full rate 
function for the probability distribution of the magnetisation, unlike earlier 
approaches (^. 28l |. Since the average over the spin variables is done over 
the Q measure, the method extends straightforwardly to arbitrary p for any 
discrete spin and for any distribution of the random helds. This is unlike 


other approaches which for p > 2, were doable only for the spherical spins |29 
Let Qi{si), for every i G N, be the probability measure on {1, —1} induced 
by the held hi + H. Then 


Q^{l) = 
Q^{-l) = 


.mi+H) 


2cosh(/3(hj + H)) ’ 

^-I3(hi+H) 

2cosh(/3(/ii -t- H)) 


( 5 ) 


We denote by Q({si}), the product measure which is a prob¬ 

ability measure on {1, —1}”. We hrst calculate an LDP with respect to 
(5({si}). We can employ Gartner Ellis Theorem for the sum, Sn = Xlili 
as it is a sum of independent random variables under the product measure 
(5({sj}). Gartner-Ellis Theorem yields that S^/N satishes LDP with respect 
to Q({sj}), with a rate function R{x) = sup^^g^ {xy — A(?/)}, where 


A(a;) = lim —hi 
N^oo N 


lim — 

V—>oo N 


N 


In 


2=1 


^x+yhi+pH _j_ ^-x-phi-pH 


2cosh(/3(hj -|- H)) 


lim {fnix) - /v(0)). 

N^oo 


( 6 ) 


Now because of the Strong Law of Large numbers, we have 


lim f^ix) ^ fix), 

N^oo 


( 7 ) 
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where /(x) is the expectation value with respect to the distribution of the 
random fields. For the bimodal distribution of the random field p{hi), given 


by 

pin = ^ [Sihi - h) + 6{hi + h )], 

(8) 

the function /(x) 

comes out to be 


fix) = 

^ln[cosh(x + (3h + /3iL)cosh(x — f3h + /3iL)]. 

(9) 

Hence we get, 

R{x) = sup {xy - k{y)] = f*{x) + /(O) 

yeR 

(10) 

where f*{x) is the Legendre-Fenchel transform of /(x). 



nx) = sup{x 2 / - fiy)}. 

(11) 


ym 


Now the probability that the sum Sn/N = x for a given /?, h and 

H, can be writtem in terms of Q({si}) as follows: 




exp(/jiVxP/p!)Q({si} ■.Y.Si = Nx) 
EL=-jvexp(/?iVxP/p!)Q({si} ; 
exp{/3NF{x))Q{{si} : 

EL=-ivexp(/?iVF(x)/2)Q({si} :J2si = Nx) 


( 12 ) 

(13) 


where, 


F(x) = 


(14) 


f3x^/p\ if |x| < 1 
f3/p\ otherwise . 

The tilted LDP ^ then yields that S^/N satisfies a LDP with respect to 
with the rate function given by 

4(x) = R{x)-F{x)-M{R{y)-F{y)} 

ySR 

= nx) + /(O) - F{x) - inf {ny) + /(O) - F( 2 /)} 


= nx) 


l3xP 

p\ 


inf < ny) 


p\ 




= sup {xy - f{y)} - —xP - inf sup {yz - f{z)} - —yP } .(15) 

yeR p'- yeR p\ 
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Let 


W{y) = xy - f{y) = xy --\og[cos\i{2{y + (3H)) + cos\i{2(3h)]. (16) 


The function W{y) attains maximum at 


y = -sinh ^ 
^ 2 


X 


1 


x^ 


ib+VT+~i^) \ - I3H, 


( 17 ) 


where a = sinh(2/Sh) and b = cosh(2/Sh). Hence the rate function in Eg. (1151) 
becomes 


jS X 

Ip{x) = log 2 — ^x^ — l3xH + — Arcsinh 


ilog 


P 

1 — x'^ 


1 — x"^ 


{b + \/l + x'^a?) 

^log (b + \/l + x‘^a?\ — inf Gi{y) (18) 
2 \ / 


where, 


Gi{y) = sup{yz-f{y)}-^yP. 
zeR P- 


At a given value of /3, h and H, the value of magnetisation m, which is the av¬ 
erage value of ■s*) averaged over all spin conhgurations, will be given 

by that value of x at which Ip{x) is minimum. Hence magnetisation is given 
by that solution of the equation Ip{x) = 0, which minimises Ip{x). Rate func¬ 
tion dehned by Eg. [THl has one local minimum at a: = 0 in the paramagnetic 
phase and develops multiple minima in the ferromagnetic phase. This rate 
function actually looks very similar to the Landau free energy functional (see 
Fig. [I]). This is not a priori obvious due to the presence of random held. 
Expanding Eg. [18] in powers of x for h = 0 and iL = 0, gives the following 
expressions for p = 2, 3 and 4 respectively: 


and 


hix) = (1 - I3)x^ -x'^ 0{x^), 

6 

hix) = x^ — —x^ + -x^ + 0{x^), 
3 6 


h{x) = x^ + 


( 2 -/ 3 ) 


12 


-h —X® -h 0(x®). 
15 


(19) 

( 20 ) 

( 21 ) 
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Figure 1: The first row shows the rate function {Ip{x)) as it changes across a continuous 
transition for p = 2. In second row we have plotted the rate function(/j,(a;)) across a first 
order transition for p = 3. 

As expected, the expansion of l 2 {x) in the absence of randomness gives 
the Landau free energy functional of the Curie Weiss model. Interestingly, 
for p > 3, even in the pure case, the /3 dependence is no longer with the lowest 
power of X. Since the expansion is not really useful to get the behaviour near 
the hrst order transition (for which we need to keep terms till all orders), 
we will not use this method to study the phase diagram of the problem 
in the presence of randomness. The full rate function contains much more 
information than just the value of magnetisation. We can easily use it to 
extract the magnetic susceptibility of the system. Magnetic susceptibility is 
the response of the system to external magnetic field (y = 

If Ip{x) gives us an equation of the form: 

( 22 ) 

where contains those terms of I'p{x) = 0, that are independent of H and 
p. Differentiating the above expression with respect to the external field H 
and then taking the limit iL —)■ 0, we get 

+ = (23) 


^ /? ^ r ,-2 • 

dx (p-2)!-^ 
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which gives. 


(24) 





We will use this expression to study the fluctuations in random fleld p-spin 
model in Section 01 

3. Phase-diagram 

For p = 2, the above rate function has already been studied and we 
will only briefly mention the results here. For p = 2 one can actually study a 
simpler function, using the fact that the global extremal points of a function 
F{x) — |||a;|p coincides with the global maximum points of |||x|p — F*{x), 
where F* is the Legendre Fenchel transform of F{x) j^. Hence the minima 
of Ip{x) will coincide with the minima of a function G{x), which is equal to: 

G{x) = — ^log(cosh[/5(x + h)]cosh[/9(x — h)]). (25) 

Interestingly this matches with the expression of free energy for random fleld 
Ising model by Salinas and Wreszinski j^. Studying the minima of this 
function, one finds that the system undergoes a continuous transition as long 
as 1 —/9sech^(/5h) < 0 and a first order transition beyond it. At equality, the 
system has a tricritical point at htp = |cosh“^y^372(~ 0.439) and {3 = 3/2. 

For p > 3, we found that the system either has first order transition to 
the ferromagnetic state or no transition at all. As shown in Fig. [H the rate 
function develops new minima discontinuously. We And that for all p > 2, 
the system fails to order into a ferromagnetic state, even at zero temperature, 
for h > 1/p!. For example, for p = 3, beyond ho ~ 0.167, the magnetisation 
in the system is always zero. Phase diagram for p = 3 and p = 4 are plotted 
in Fig. |2j 

4. Magnetic susceptibility 

We want to now look at the fluctuations in the magetisation, as given by 
the magnetic susceptibility of the system. Note that x(a;,/?, h) is a function 
of magnetisation itself, and hence one needs to get the value of x for which 
Ip{x) is minimum, and then substitute it in the equation of x(x, h, (3) to get 
the physically acceptable value of magnetic susceptibility at a given h and (3. 

4-1. p=2 

For p = 2, in the paramagnetic phase, x(0, /3, h) has a particularly simple 
expression: 
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Figure 2: Phase diagram for p = 3(solid line) and p = 4(dotted line). Below the line system 
is in ferromagnetic phase and above it, it is in paramagnetic phase. For h > ho ~ 0.167 
for p = 3, and for h>ho~ 0.041 for p = 4, the model has no ferromagnetic phase, i.e, it 
has zero magnetisation for all values of /3 for h > ho 


x{0,P,h) = 


2/9 


(26) 


1 — 2/9 + cosh(2/9/i) 

For h = 0, the above equation gives a divergence in x at /9c = 1. But for 
h 7 ^ 0, it gives two values of (3 where x(0, /9, h) diverges for all h < 0.447. We 
hnd that for h < htp, the lower value of (3 is the physically relevant one, as 
beyond that value of f3, magnetisation is no longer zero. For htp < h < 0.447, 
the (3 at which the transition occurs is less than the values of [3 that one 
obtains by equating x(0,/3, h)~^ to zero. 

For p = 2, near the continuous transition (for h < htp ^ 0.439), we can 
study the magnetic susceptibility by keeping terms only till in Eq.(|22]) by 
expanding J-{x, f3, h) around a; = 0. We get 

3 

—/377+—(1 — 2/9+ cosh(2/9h)) — —[cosh^(/9h){—2 + cosh(2/9h)}] = 0. (27) 
2 3 

Differentiating with respect to H and taking the limit 7? ^ 0, we get: 


x(0,/3,h) = 


2/9 


[1 — 2/9 + cosh(2/9h)] — 2a:^cosh"^(/9h)[cosh(2/9h) — 2] 


(28) 


Putting h = 0 in the above equation, one gets y ~ for /9 < 1 and 
X ~ for /9 > 1, as expected. Hence, in the absence of random held, we 

reproduce the results of standard mean-held theory. 

For small values of h, we can further expand the cosh term and in that 
regime, the magnetisation is given by: 


m = 0 


, for 13 < (3c 
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V3(-l + /g-/g2/.2) , for /3 > /3„ 


(29) 


where 

Substituting in the expression of x, we get, for (3 < (3c 


and for (3 > (3c 


X{0,b, h) ~ 


l-(3 + (3^h^' 


X{m,b,h) ~ 


2(l-/3 + /?2h2)' 


(30) 


(31) 




Figure 3: Plot of x as a function of /3 for h > htp 

For h > htp{^ 0.439), the system undergoes a hrst order transition. For 
h > htp, since the system has a hrst order transition we have to solve Eq. |22] 
and Eq. [2l] numerically to obtain the magnetic susceptibility. We have plot¬ 
ted magnetic susceptibility and magnetisation for a few representative values 
of h in Fig 3. We hnd that for htp < h < 0.447, y is rapidly increasing as one 
approaches the transition point from the high temperature phase. On the 
other hand, for h > 0.447, y peaks at a value oi (3 < (3^ and then decreases 
as it approaches (3c- This difference in behaviour of magnetic susceptibility 
arises due to the fact that for htp < h < 0.447, even though the system 
has a hrst order transition, the function x(0,/d,h) still shows divergence at 
two values of (3. In this regime, none of them is relevant as none of them 
correspond to the transition point and occur well inside the ordered regime. 
But due to their presence, for htp < h < 0.447, y still increases rapidly as 
one approaches (3c- 

4-2- p > 2 

For p = 3 and p = 4, we numerically evaluate the expression in Eq. 
1241 and hnd that the x(0,/3, h) is a convergent function for all values of (3 
and h- Magnetic susceptibility shows a discontinuity, as expected, at the 
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transition point. In Figs. 4 and 5, we have plotted magnetic snsceptibility 
and magnetisation for a few representative valnes of h for p = 3, and p = 
4, respectively. For low /?, below the transition, y oc /9 as expected for 
a paramagnet. Interestingly for h > 1/p!, even thongh the system does 
not order ferromagnetically, magnetic snsceptibility shows a non monotonic 
behavionr as a fnnction of (3 as shown in Fig HI Similar behavionr of magnetic 
snsceptibility is also seen for p = 4 beyond h > 0.041 (Fig. H]). 




Figure 4: m and y for p = 3 for h < 0.167 and for h > 0.167. 




Figure 5: m and y for p = 4 for < 0.041 and for h > 0.041. 


5. Conclusion 

The nature of phase transition of p-spin random held model is qualita¬ 
tively different for p = 2 and for p > 2. For p = 3, the random held model 
has been solved for spherical spins earlier j2^. We solve it here for Ising 
spins, and hnd that for all p > 2, for high enough strength of the random 
held, Ising spins do not order into a ferromagnetic phase, even at zero tem¬ 
perature. Magnetic susceptibility for these strengths of helds also shows a 
non monotonic behaviour, peaking at a hnite f3. The method can be used 
starightforwardly to study other random distributions, though for p > 3, we 
do not expect any qualitative diherence in the behaviour of the model for 
diherent distributions. 
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We also notice that the nature of fluctuations near the transition point 
is different for the hrst order transition when p = 2, from that of p > 3. For 
p = 2 there exists a range of h for which even though the system undergoes 
a hrst order transition, x increases rapidly as one approaches the transition. 
It would be interesting to see if the similar behaviour continues in hnite 
dimensions. In three dimension, RFIM is known to have strong hnite size 
ehects which makes it difficult to acertain the nature of transition for the 
higher strengths of the random held jsl, 1^. Besides that, Eq. [18] can also give 
us some indications for writing down the correct held theory for the model. 
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